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We study four different models of Chern insulators in the presence of strong electronic repulsion 
at partial fillings. We observe that all cases exhibit a Laughlin-like phase at filling fraction 1/3. We 
provide evidence of such a strongly correlated topological phase by studying both the energy and the 
entanglement spectra. In order to identify the key ingredients of the emergence of Laughlin physics 
in these systems, we show how they are affected when tuning the band structure. We also address 
the question of the relevance of the Berry curvature flatness in this problem. Using three-body 
interactions, we show that some models can also host a topological phase reminiscent of the v = 1/2 
Pfaffian Moore-Read state. Additionally, we identify the structures indicating cluster correlations 
in the entanglement spectra. 
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I. INTRODUCTION 

First introduced by Haldane in his 1988 paper^, the 
Chern insulator is the simplest example of a topologi- 
cal insulator. The topological insulating phase is charac- 
terized by the (non-zero) Chern number of the occupied 
band, and it exhibits a non-zero integer Hall conductance 
a xy at zero magnetic field. The Hall conductance can be 
attributed to the metallic edge states protected by the 
non-trivial topology of the bulk. Interest in such states 
was boosted by the theoretica l pr edictiorPHU and subse- 
quent experimental realization^!! of time-reversal topo- 
logical insulators. In its simplest form in two dimen- 
sions, this type of insulator consists of two copies (spin up 
and down) of a Cher n in sulator. Comprehensive studies 
have been carried ouiP^, and a complete classification of 
topological band theories with time-reversal and charge- 
conjugation sy mmetries in all dimensions has been es- 
tablishecP^HH]. More than two decades after its first in- 
troduction, however, the theoretical study of topologi- 
cal insulators is still mostly limited to the single-particle 
regime, with interaction effects playing a subleading role. 

Recent developments have shown that strong electronic 
interactions open up interesting new possibilities. Novel 
phases resembling the fractional quantum Hall (FQH) 
effect have been identified in lattice models at zero mag- 
netic field. Several authors have reported the discovery of 
the 'fractional Chern insulator' (FCI)PHlI]_ They found a 
FQH phase of interacting electrons on a checkerboard lat- 
tice at filling 1/3 and zero magnetic field. The feature of 
this phase is an almost 3-fold degenerate incompressible 
ground state with Hall conductance a xy = 1/3, reminis- 
cent of the FQH Laughlin state on a torus. A similar 
state has been found on a triangular lattice as welP^. 
In par allel, F CI phases of interacting bosons have been 
realizecf^E^, a time-reversal-symmetric fractional topo- 
logical hquid state has been constructed from two copies 
of FCI 19 * 2 ° I, and an integer quantum Hall effect has been 
found in a half-filled Hubbard model in coexistence with 
an Ising ferromagnetic ordeJ23. 



The resemblance of the FCI phase to the conventional 
FQH effect has been justified from several perspectives. 
First, the FCI ph ase ha s quasiparticle excitations with 
fractional statistic d 15 * 22 l, and the quasiparticles are gov- 
erned by the admissible count ing rules first found in the 
study of the FQH effectP^. Second, in the limit of 
long wavelength and uniform Berry curvature, the pro- 
jected single-particle density operators form a closed Lie 
algebr£p322]. This algebra has the same structure as the 
Girvin-MacDonald-Platzman algebra of magnetic trans- 
lations and projected density operators in the FQH ef- 
feclpS Third, Wannier functions maximally localized in 
one dimension have been explicitly constructed from the 
single-particle states of a Chern bancP^. Finally, var- 
ious constructions of FQH-analogue wave-functions for 
the ground state of FCI have been proposed, using either 
Wannier functions^ or a parton approaclP^HH] 

In this paper we show that the FCI phase is present 
in the Haldane model on the honeycomb lattice^, in 
a two-orbital model that resembles half (spin-up) of 
the Mercury- Telluride two-dimensional topological insu- 
latoiP, in the Kagome lattice model with spin-orbit cou- 
pling^!, anc i i n the spin-polarized ruby lattice modeP^. 
These models allow us to study FCI in different physi- 
cal situations such as different lattices or different num- 
ber of Bloch bands. Working in a single flattened band, 
we find in each model a more or less robust Laughlin 
FQH phase at 1 /3 filling in the presence of repulsive two- 
body nearest-neighbor interactions. The 3-fold degener- 
ate ground states are separated from the excited states by 
a finite gap and flow into each other upon flux insertion 
with a period of 3 fluxes. This signals a Hall conductance 
a xy — 1/3. We identify hallmarks of fractional excita- 
tions of the Laughlin 1 /3 universality class in the energy 
spectrum as well as the entanglement spectrum. We then 
discuss the stability of the topological ground state un- 
der parameter variation, and test its correlation with the 
anisotropy of the Berry curvature. We also show that 
another FCI phase reminiscent of the Pfaffian Moore- 
Read FQH state^ is present in the Kagome and ruby 
lattice models at half filling of the valence band. Finally, 
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FIG. 1: The Haldane model on the honeycomb lattice. A and 
B sublattices are colored in red and green, respectively. The 
lattice translation vectors are bi = a.2 — a$ and b2 = a-s — a± . 
The hopping between NN is t\ , and the hopping between NNN 
is t^e 1 in the direction of arrows. The density-density repul- 
sion between nearest neighbors is depicted in gray ellipses. 



we highlight several structures in the higher levels of the 
particle entanglement spectru m 1 35 * 36 ! of the ground state 
at filling 1/3 that may serve as a hint for the stable ex- 
istence of other FQH states at other fillings, such as the 
Read-Rezayi series^. 



II. HALDANE MODEL 

The Haldane modef^ is the first studied example of a 
topological insulator. We would like to see if this model 
can host fermionic FCI phases (the bosonic version has 
recently been reported in Refs. ll7lT5)l . We adopt the hon- 
eycomb lattice layout from Ref. [TBJ As shown in Fig. [T] 
the two sublattices A and B are connected by the vec- 
tors ai = (0, -1), a 2 = (V3/2, 1/2), a 3 = (-s/3/2, 1/2). 
We define the lattice translation vectors bi = a 2 a 3 , 
b 2 = a 3 — ai. The Haldane modeP has real hopping 
amplitude t\ between nearest neighbors (NN), complex 
hopping amplitude t 2 e ±l< ^ between next-nearest neigh- 
bors (NNN) , and an inversion-breaking sublattice poten- 
tial M. 

After a Fourier transform to the first Brillouin 
zone and a gauge transform on the B sublat- 
tice Vk.s -> e ik '( bl+2b2 )/ 3 i/'k„B, the single-particle 
Hamiltonian can be put in Bloch form as H = 

Sk(Wc,A'V , I, B )MKKV'k,A,'i/'k,B) T - Here the lattice mo- 
mentum k = (k • bi,k • b 2 ) = (k x ,k y ) is summed over 
the first Brillouin zone, and the /i(k) matrix can be ex- 
pressed in terms of the identity and 3 Pauli matrices, 



h(k) = d I + J2i di<Ji, where 

do = 2i 2 cos 0[cos k x + cos k y + cos(k x + k y )] , (1) 

d% = £i [1 + cos(k x + k y ) + cos k y ], 

d y = ti[— sm(k x + k y ) — sinfcy], 

d z = M + Iti sin0[sinfc. r + sink y — sm(k x + ky)]. 

The single-particle Hamiltonian has inversion symme- 
try at M = and the 3-fold point group symmetry of 
the honeycomb lattice. At AI = 0, inversion exchanges 
the two sublattices and transforms the annihilation op- 
erators by (V'k,yi,'0k,B) T -> o- a .(-0_ kiA ,'0-k,s) T - The 3- 
fold rotation generates the cyclic permutation of the lat- 
tice translation vectors bi — > b 2 > — bj b 2 > hi on 
each sublattice and thus transforms the wave vectors by 
(k x ,k y ) — > (k y ,—k x — k y ). Therefore, the Bloch Hamil- 
tonian has the following two symmetries: 

h(k x ,k y ) = cr x h(-k x , -k y )<j x , (2) 
h(k X: ky) = U {k x ,ky)h(ky, k x ky)U {k x , ky) , 

where U(k x , k y ) is a diagonal 2x2 unitary matrix, with 
[1, e l ( k *+ k v)] on the diagonal. When the system is put on 
the lattice of finite size N x x N y with periodic boundary, 
the 3-fold rotation symmetry is lifted, unless N x — N y . 

To focus on the effect of interactions without being 
distracted by single-particle dispersion, we always take 
the flat-band limit of the insulator, i.e. replace the 
original Bloch Hamiltonian h(k) — P„(k)P„(k) by 
^] n £ n (0)P n (k), where P„(k) is the projector onto the 
n-th band. We then send the band gap to infinity and 
work directly in the lowest band, in the same spirit of the 
lowest Landau level projection routinely adopted in the 
FQH literature. We then add a density-density repulsion 
between nearest neighbors. Since a flattened band does 
not provide an energy scale, we are free to set the inter- 
action strength to unity. After the aforementioned gauge 
transform on ipKB, the interaction term can be written 
in the sublattice basis as 

i E <A< B V'k 2 B^k 1 Ac 1 o + d k:-k 3 -k 4 ^ 1 k 2 k3k 4 , (3) 

where 

Vk lk2k 3k 4 = 1 + e l ^~ k ^ + e i(k2-k 4 )(b 1 +b 2 ) ! (4) 
as illustrated in Fig. [T] 

A. Ground state at 1/3 filling 

We diagonalize the interacting Hamiltonian in the flat- 
tened lowest band at filling 1/3. We show the energy 
spectrum of N = 8, 10, 12 particles on the N x x N y = 6 x 
Y lattice in Fig. 2 The calculations are performed with 
(ti, i 2 , M, 4>) = (171, 0, 0.13). The particular choice of pa- 
rameters will be discussed later. In the three cases (N — 
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FIG. 2: Low energy spectrum of the Haldane model with 
N — 8 (marked by red crosses), A^ = 10 (green plus signs), 
and N — 12 (blue triangles) particles on the N x x N y — 6 x ^■ 
lattice, with energies shifted by Eq, the lowest energy for each 
system size. We only show the lowest excited level in each 
momentum sector in addition to the 3-fold ground state. 
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FIG. 3: Energy gap AE for different system sizes and aspect 
ratios of the Haldane model. The gap is defined as the energy 
difference between the first excited state and the highest of the 
3-fold ground states. A more precise definition of AE is given 
II D| In each case, N x = 3N/N y . The negative 



in Section 

values of AE in the N y = 3 group mean that the actual 
ground state is not in the momentum sectors predicted by 
the theoretical counting rule, i.e. the gap above the supposed 
topological ground state closes. 



8, 10, 12), a 3-fold degenerate ground state is seen at total 
momenta {(0, 0), (2, 0), (4, 0)}, {(1, 0), (3, 0), (5, 0)}, and 
{(0, 0), (0, 0), (0, 0)}, respectively. This agrees perfectly 
with the (1, 3)-admissible counting proposed in Ref. [Tjjl 
and recently developed in Ref. [251 The finite-size scaling 
of the energy gap AE is shown in Fig. [3] The energy gap 
does not seem to remain open in the thermodynamic limit 
even if the aspect ratio N x /N y remains finite. Nonethe- 
less, as discussed in the following sections, we find solid 
evidence for the topological nature of the ground state; 
the detailed investigation into the energetics will be pre- 
sented in future work. 



FIG. 4: Spectral flow of the low-lying states of the Haldane 
model with A^ = 10 particles on the N x x N y = 6x5 lat- 
tice upon flux insertion along the x direction. ~y x counts the 
number of fluxes inserted. The 3-fold ground states flow into 
each other, and do not mix with higher states during flux in- 
sertion. It takes 3 full fluxes for the 3-fold states to return to 
the original configuration (inset). 



The three degenerate ground states exhibit spectral 
flow upon flux insertion. For example, inserting a unit 
flux in the x direction shifts the single particle momen- 
tum k x — > k x + 2ir~/ x , with ~f x going from to 1; this 
induces the spectral flow within the 3-fold ground state 
(Fig. [1]). Upon insertion of 3 full fluxes, the 3 degener- 
ate states restore their original configuration. Given the 
unit Chern number of the valence band, we conclude the 
system has Hall conductance a xy = 1/3. 

We have also checked the effect of density-density re- 
pulsion between NNN on the 1/3 phase. Overall, this 
additional interaction term weakens the Laughlin-like 
phase. These results are in agreement with those of 
Sheng et alP^ for the checkerboard model. 



B. Quasihole excitations 

Three-fold degeneracy alone is insufficient to fully es- 
tablish the observed ground state as a FQH state; it could 
very well be an imprint of a lattice charge density wave. 
To rule out this alternative we need to study the exci- 
tations of the system. We decrease the filling v from 
1/3 and check for quasiholes in the energy spectrum. In 
Fig. [5] we show the energy spectrum of N = 7 particles 
on the N x x N y = 6x4 lattice, and in Fig. [6] the energy 
spectrum of N = 9 particles on the N x x N y = 6x5 
lattice. (We run into convergence problems at A^ = 11.) 
These configurations have the same lattice geometry as 
the corresponding ground states shown earlier, but with 
one electron removed. An energy gap is clearly visible 
in the spectrum, and the low-energy part has the same 
counting in each momentum sector as predicted by the 
(1, 3)-admissible rule. This further substantiates that the 
ground state observed at filling 1/3 indeed has the basic 
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FIG. 5: Low energy spectrum of the Haldane model with 
TV = 7 particles on the N x x N y =6x4 lattice. The number 
of states below the gray line is 12 in each momentum sector, 
in agreement with the (1, 3)-admissible counting rule. 




FIG. 6: Low energy spectrum of the Haldane model with 
TV = 9 particles on the TV,, x N y =6x5 lattice. The number 
of states below the gray line is 19 in momentum sectors with 
k x — 0,3 and 18 elsewhere, in agreement with the (1,3)- 
admissible counting rule. 



features of the Laughlin FQH stat* 
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C. Entanglement spectrum 



Recent development d 35 * 3 ^ * 3 ^ showed that the excita- 
tions in FQH systems are manifested in the entangle- 
ment between particle j 35 l 40 l Ml i n the many-body ground 
state. Using this alternative probe, we provide further 
evidence that the ground state at filling 1/3 is a FQH 
Laughlin state. This tool is highly valuable in the present 
case since no overlap with model wave functions can be 
computed: despite several proposal^SHaiJ concrete ex- 
pressions for the model wave functions have not been 
established for FCI. 

Specifically, we cut the system in the way described in 
Ref. |3H and further used to look at the FCIs in Ref. IT51 



FIG. 7: Particle entanglement spectrum of the ground state 
of the Haldane model of TV = 12 particles on the N x x N y = 
6x6 lattice, with TVs = 8 particles traced out. The number of 
states below the gray line is 741 in momentum sectors where 
both k x and k y are even and 728 elsewhere, in agreement with 
the (1, 3)-admissible counting rule. The system size is exactly 
the same as Fig. 13 of Ref. 115) 



We divide TV particles into two subsystem of TV^ and TV# 
particles, and trace out the degrees of freedom carried by 
the Nb particles. The eigenvalues e - '' of the resulting re- 
duced density matrix pa — TtbP — T^B\ip)(ip\ define the 
particle entanglement energies £. For degenerate ground 
states, we form the density matrix as an incoherent sum 
with equal weights^ P = \ J2i IV'iXV'il- Then, the en- 
tanglement energy levels £ can be displayed in groups 
marked by the total momentum {k x , k y ) of the TV^ par- 
ticles. A typical case is shown in Fig. [7] The spectrum is 
very similar to what was found in the checkerboard lat- 
tice modeP=. We observe a clear, although narrower, en- 
tanglement gap in the spectrum; the counting of the en- 
tanglement energy levels below the gap matches (in each 
momentum sector) the (1, 3)-admissible quasihole count- 
ing^ of TV^ particles on the TV^ x TV y reciprocal lattice. 
We have checked all values of Na manageable by current 
computers and have found perfect agreement with the 
counting principle in all cases. Given the vast difference 
between the checkerboard lattice model and the Haldane 
model (lattice symmetry, coordination number, flux dis- 
tribution, etc.), the similarity in the entanglement spec- 
trum is surprising. Compared with the relatively fuzzy 
quasihole energy spectrum, the ground state entangle- 
ment spectrum turns out to be a more reliable alternative 
route to probe the physics of fractional excitations. 



D. Parameter dependence 

As analyzed by Haldane 1 , strong inversion breaking 
eventually overcomes the non-trivial topology; the Chcrn 
number vanishes when \M\ > 3\/3|£2 sin0|. We now 
study whether strong interactions would further shrink 
the volume of the topologically non-trivial phase in pa- 



FIG. 8: The distribution of the energy gap AE in the 
(4>, M) plane for the Haldane model of N = 8 particles on 
the N x x N y =6x4 lattice. The color code for AE is shown 
by the bar on the right, and the region with the gap closed 
is marked in gray. The bold red line M = 3v3sin0 sepa- 
rates the topologically non-trivial sector of the single-particle 
problem from the trivial one. 



rameter space. The density-density interaction projected 
to the flattened valence band has two parametric degrees 
of freedom (M/ti,t2sixi<f>/ti), which effectively change 
the Berry curvature in the model. Without loss of gener- 
ality, we fix ti = t% = 1 and study the region with M > 
and <p G [0, 7r/2]. There is no clear boundary between dif- 
ferent interacting phases in finite-size numerics. We thus 
need a quantitative characterization of the similarity to 
the ideal FCI state. 

We denote the (1, 3)-admissible counting of N parti- 
cles on the N x x N y lattice in momentum sector (k x , k y ) 

by n k ^ k °" Nv . We denote by A the collection of the low- 
est n k j k ' Ny states in each momentum sector, and by 

A the collection of all the other states. If the system is 
in a well-developed FCI state, the collection A is the 3- 
fold degenerate ground state, while the collection A con- 
tains the excited states. The energy gap is thus AE — 
mhiE^ — max.E_4, and the energy spread of the ground- 
state manifold can be defined as SE = max Ea — min Ea . 
Further, we calculate the entanglement spectra of the de- 
generate ground state for various partitions (Na, Nb) of 
N particles. For each Na, we denote by A the collec- 
tion of the lowest n k A ' k °" Ny entanglement energy levels 

in each momentum sector, and by A all the other entan- 
glement energy levels. We define the entanglement gap 
A£ = min £^ — max£&- The parameter sets with large 
AE, small SE and large A£ are likely to host a FCI state. 

In Figs, m |9j [lOj we show the distribution of AE, 5E, 
A£ over the (4>, M) plane. Combining the three plots, we 
find that the region enclosed by \M\ < 3\/3|*2 sin 0| and 
\M\ < 3\/3|*2|(0.4 - I sin 0|) has large AE and A£, and 
small SE, and hence is likely to support a FCI phase. 



FIG. 9: The distribution of the energy spread SE in the 
((f), M) plane for the Haldane model of N = 8 particles on the 
N x x N y = 6x4 lattice. The color code for SE is shown by 
the bar on the right. The bold red line M = 3\/3sin0 sepa- 
rates the topologically non-trivial sector of the single-particle 
problem from the trivial one. 




particles traced out. The color code for A£ is shown by the 
bar on the right, and the region with the gap closed is marked 
in gray. The average spacing of the entanglement levels below 
the gap is at most 0.05 (not shown in the figure); hence the 
entanglement gap should be considered widely open as long as 
A£ 2> 0.05. The line with A£ = 0.5 is shown in magenta. The 
bold red line M — 3\/3sin</!> separates the topologically non- 
trivial sector of the single-particle problem from the trivial 
one. 



E. Berry curvature variation 

A few authors pointed out that a Chern band sh ares 
an important feature with the Landau leveP322E6|_ j n 
the limit of long wave-length and uniform Berry curva- 
ture, the projected density operators form a closed Lie 
algebra. This algebra has the same structure as the 
Girvin-MacDonald-Platzman algebra of magnetic trans- 
lations^, with the Berry curvature taking the role of the 
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FIG. 11: The distribution of the standard deviation of the 
Berry curvature erg in the (<j>, M) plane for the Haldane model 
of N = 8 particles on the N x x N y = 6x4 lattice. The bar on 
the right shows the color code for the value of 2ir ae- The bold 
red line M = 3\/3sin<^ separates the topologically non-trivial 
sector of the single-particle problem from the trivial one. We 
focus on the non-trivial sector. The strong fluctuation of the 
Berry curvature very close to the boundary of the non-trivial 
sector is not shown. 



uniform magnetic field. They thus argued that the de- 
velopment of the FCI phase is driven by this algebraic 
structure. This picture suggests a negative correlation 
between Berry curvature fluctuations and the propensity 
towards a FCI phase. This is indeed observed in finite- 
size numerics, as we describe below. 

We measure Berry curvature fluctuations by the sim- 
plest possible option, its standard deviation trg in units 
of the average Berry curvature \B\ = 1/2tt. The distri- 
bution of 2it<jb is shown in Fig. |11| Comparing this with 
the patterns in Figs. [81 [9 10 we find a clear correlation 



between erg and the three measures of the propensity 
towards the FCI phase over the full range of the param- 
eter scan. In particular, the optimal parameter region 
supporting a robust FCI phase coincides with the part 
of the parameter space with least fluctuating Berry cur- 
vature. While this agrees with the picture we expect 
from the algebraic structur e! 22 * 25 * 26 * 4 ^, we stress that the 
fluctuation at the optimal parameters is still quite sig- 
nificant, with its standard deviation comparable to the 
mean value. The FCI is apparently more robust to Berry 
curvature fluctuations than expected. This is crucial for 
two-band Chern insulators since their Berry curvature 
cannot be completely flat due to the no-hair theorem^. 



III. TWO-ORBITAL MODEL 

We now turn to a simpler model with two orbitals on 
each site of a square lattice. This model represents the 
spin-up half of the Mercury- Telluride tw o-d imensional 
topological insulatoJ^. As shown in Fig. 12 there are 
two orbitals A, B on each site, with energy difference 
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FIG. 12: The two-orbital model on square lattice. The two 
orbitals A,B are colored in red and green, respectively. The 
intra-orbital NN hopping amplitude is t2 between the A or- 
bitals, and —t2 between the B orbitals. The inter-orbital NN 
hopping amplitude is ±it\ in the x direction along the arrows, 
and ±ti in the y direction. The on-site Hubbard repulsion is 
depicted by the ellipse in gray. 



2M . The intra-orbital NN hopping has amplitude ±t2, 
while the inter-orbital NN hopping has amplitude ±i£i 
and ±ii in the x and y directions, respectively. These 
amplitudes have been specifically designed such that the 
single-particle Bloch Hamiltonian takes the form H = 
Ek(V'I,A'' ! /'k,B) /l ( k )(V'k^,'0k,s) T , where 

h(k) — 2ti sin k x a x + 2t\ sin k y a y 

+ [M — 2t2(cosk x + cosk y )](x z . (5) 

Similar to the case of the Haldane model, we can flatten 
the Bloch bands using projectors. We then add Hubbard 



inter-orbital repulsion to each site, as shown in Fig. 
In momentum space, the interaction term reads 
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N ^ 



°ki +k 2 -k 3 -k 4 Vk 3 A ^k 4 



jV'kaS^: 



(6) 



Even though at single-particle level a two-orbital-per-site 
model is equivalent to a model with two sites in each 
unit cell, this situation changes when interactions are in- 
cluded. Notice that the interaction here (Eq. [6]) has a 
different form factor than the one in Eq. [4j 

We diagonalize the interacting Hamiltonian in the flat- 
tened lowest band at filling 1/3 and fe/ti, M/ti) — 
(1,2). We show the energy spectrum of N = 8. 10 par- 

(We 
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tides on the N x x N y = 6 x | lattice in Fig. 
run into serious convergence problems at N = 12.) In 
the two cases (N = 8,10), a 3-fold degenerate ground 
state is seen at total momenta {(0,0), (2,0), (4,0)} and 
{(1, 0), (3, 0), (5, 0)}, respectively. This agrees perfectly 
with the (1, 3)-admissible counting proposed in Ref. IT51 
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FIG. 13: Low energy spectrum of two-orbital model for N = 
8 (marked by red crosses) and N = 10 (green plus signs) 
particles on (N x , N y ) = (6, N/2) lattice, with energies shifted 
by Eo, the lowest energy for each system size. We only show 
the lowest excited level in each momentum sector in addition 
to the 3-fold ground state. 
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FIG. 14: Energy gap AE of the two-orbital model for differ- 
ent system sizes and aspect ratios. In each case, N x = 3N/N y . 
The negative values of AE in the N y = 3 group mean that the 
actual ground state is not in the momentum sectors predicted 
by the theoretical counting rule. 



l22l As shown in Fig. [Mj the energy gap AE remains 
open and scales to a finite value in the limit of N — > oo 
with N x /N y finite. The three degenerate ground states 
exhibit spectral flow upon flux insertion with a period 
of 3 fluxes (see Fig. 15 1. This shows that the system 



has Hall conductance a xy = 1/3. We study the quasi- 
hole excitations through ground state entanglement. In 
Fig. [16] we observe an entanglement gap in the spec- 
trum. The counting of the entanglement energy levels 
below the gap again matches in each momentum sector 
the (1, 3)-admissible counting. This shows that the exci- 
tations in the ground state of the two-orbital model have 
the same counting as that of Abelian fractional statistics 
1/3 quasiholes. We thus conclude that the ground state 
is a FQH Laughlin state. 



FIG. 15: Spectral flow of the low-lying states of the two- 
orbital model with N = 10 particles on the N x x N y = 6x5 
lattice upon flux insertion along the x direction. y x counts 
the number of fluxes inserted. The 3-fold ground states flow 
into each other, and do not mix with higher states during flux 
insertion. It takes 3 full fluxes for the 3-fold states to return 
to the original configuration (inset). 
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FIG. 16: Particle entanglement spectrum of the ground state 
of the two-orbital model of N = 10 particles on the N x x N y = 
6x5 lattice, with Nb = 5 particles traced out. The number 
of states below the gray line is 776 in momentum sectors with 
k y = and 775 in all the other sectors, in agreement with the 



analytical counting rule. 



The two-orbital model has topologically non-trivial 
bands when \M\ < 4 j ^2 ] - We now move the system away 
from the point (t 2 ,M) = (1,2) and probe this parameter 
region. In Figs. [T7j [THJ [l9j we display the distribution of 
the energy gap AE, the energy spread 8E, and the en- 
tanglement gap A£ of the ground state. Compared with 
the Haldane model, the situation is more complicated. 
The maximum of the energy gap AE does not coincide 
with the minimum of the energy spread 5E. Rather, 
the region with large gap AE tends to have large spread 
SE as well. The peak of the entanglement gap A£ is 
close neither to the maximum of AE nor to the min- 
imum of SE. Varying (M, t2) scans over the manifold 
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FIG. 17: The distribution of the energy gap AE in the 
(t2,M) plane for the two-orbital model of N = 8 particles 
on the N x X N y = 6 X 4 lattice. The color code for AE is 
shown by the bar on the right, and the region with the gap 
closed is marked in gray. The bold red line M = 4t2 sepa- 
rates the topologically non-trivial sector of the single-particle 
problem from the trivial one. 
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FIG. 18: The distribution of the energy spread 8E in the 
(t,2, M) plane for the two-orbital model of N = 8 particles on 
the N x x N y = 6x4 lattice. The color code for SE is shown by 
the bar on the right. The bold red line M = 4t2 separates the 
topologically non-trivial sector of the single-particle problem 
from the trivial one. 



of the two-orbital interacting Hamiltonians projected to 
the lowest band. The lack of correlation between AE, 
SE and A£ suggests that the distance from this manifold 
of Hamiltonians to the Laughlin model Hamiltonian is 
quite large. Further, we check the correlation between 
Berry curvature fluctuations and the propensity towards 
the FCI phase. The variation of the curvature fluctuation 
is shown in Fig. [20] We observe weak correlation between 
erg and SE and find little association between erg and AE 
or A£. It is possible that erg is not a very good measure 
of the Berry curvature fluctuation. We note that in this 
model the interaction before projection is constant. We 
leave in-depth study of this issue for future work. 



FIG. 19: The distribution of the entanglement gap A£ in the 
(t2, M) plane for the ground state of the two-orbital model of 
N = 8 particles on the N x X N y = 6 X 4 lattice, with Nb = 5 
particles traced out. The color code for A£ is shown by the 
bar on the right, and the region with the gap closed is marked 
in gray. The average spacing of the entanglement levels below 
the gap is at most 0.05 (not shown in the figure); hence the 
entanglement gap should be considered widely open as long 
as A£ ^> 0.05. The line with A£ = 0.5 is shown in magenta. 
The bold red line M = 4t2 separates the topologically non- 
trivial sector of the single-particle problem from the trivial 
one. 




FIG. 20: The distribution of the standard deviation of the 
Berry curvature ob in the (i,2,M) plane for the two-orbital 
model of N = 8 particles on the N x x N y = 6x4 lattice. The 
bar on the right shows the color code for the value of 2na8- 
The bold red line M — 4t2 separates the topologically non- 
trivial sector of the single-particle problem from the trivial 
one. We focus on the non-trivial sector. The strong fluctua- 
tion of the Berry curvature very close to the boundary of the 
non-trivial sector is not shown. 



IV. KAGOME LATTICE MODEL 

We study the Kagome lattice model built by Tang 
et alP^l. As shown in Fig. 21 the lattice is spanned 
by the translation vectors b! and b 2 , and it consists 
of three sublattices A,B,C. In Ref. [3H the single- 
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FIG. 21: The Kagome lattice model with NN spin-orbit 
coupling. The three sublattices A,B,C are colored in red, 
green, blue, respectively. The lattice translation vectors are 
bi = 2ai, Y>2 = 2a2. The hopping amplitude between NN 
is ti — i\i in the direction of arrows. The two-body and the 
three-body interactions between NN are illustrated by the 
gray ellipses and triangles, respectively. 



particle model has been studied both without and with 
NNN hopping terms. We have looked at the effect of 
interactions in both cases. But for sake of simplic- 
ity we focus mostly on the case without NNN hop- 
pings. The hopping amplitude between NN is t\ ± 
After a Fourier transform and a gauge transform 
V>k,B -> <Ak,Be- lk ' bl / 2 , Vk,c -> <Ak,ce~ lk ' b2/2 , the single- 
particle Hamiltonian can be cast in Bloch form as H = 
Ek(V , k,A'^k,B)V , k,c)M k )(V'k,A, V>k„B,V'k,c) T - Here the 
lattice momentum k = (kb l5 k-b2) = {k x , k y ) is summed 
over the first Brillouin zone, and the h(k) matrix reads 



h(k) = -ti 



1 + e 



-ifc, 



1 



e 



- e 




1 



1 + e i ( k: *~ ky } 

e i(k y -k m ) o 



1 



-ik x 



e -"*x _l_g- 

— 1 — e ifex 1 + e i ^ kx ~ ky ^ 

I _|_ e ik y —I e i(ky-K) 



(7) 



The three Bloch bands can be flattened using the projec- 
tor method detailed for the Haldane model in Section HT1 
We focus on the lowest band. Unless specified otherwise, 
the numerical calculations shown below are performed at 
Ai = t\ as discussed in the original papepS The lowest 
band has unit Chern number. 
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FIG. 22: Low energy spectrum of the Kagome lattice model 
with N = 8 (marked by red crosses), iV = 10 (green plus 
signs), and iV = 12 (blue triangles) particles on the N x x 
N y — 6 x N/2 lattice, with energies shifted by Eq, the lowest 
energy for each system size. Also shown (magenta squares) 
is the spectrum of N = 12 particles on the N x x N y =6x6 
lattice in the alternative Kagome model with NNN hoppings 
at the parameters suggested in the original paper^J namely 
(Ai,A 2 ,t 2 ) = (0.28,0.2,-0.3). 



A. Filling 1/3 

We fill the flattened lowest band to filling 1/3, and add 
density-density repulsion between NN. After the gauge 
transform, the interaction term reads 



A,B,C 

a</3 



k 2 k 3 k 4 ' 



(8) 

where the sublattice indices (a, /3) are summed over 
(A, B), (B,C), (C, A), and the interaction factors are 



v ab _ 1 

^kjkakaki - 1 

yBC _ i 

''kjkakaki - 1 

yCA _ i 

^kjkakaki - 1 



e 



-i(k 2 -k 4 )-bi 
i(k 2 -kt).(bi-b 2 ) 



(9) 



3 i(k 2 -k 4 )-b 2 



The 6 terms are illustrated by the 6 ellipses in Fig. [21] 
We show the energy spectrum of = 8, 1(^12 parti- 



cles on the N x x N y — 6 x y lattice in Fig. 22 In the 
three cases, a 3-fold degenerate ground state is seen at 
total momenta {(0, 0), (2, 0), (4, 0)}, {(1, 0), (3, 0), (5, 0)}, 
{(0, 0), (0, 0), (0, 0)}, respectively. Again, this agrees per- 
fectly with the (1, 3)-admissible counting proposed in 
Ref. 1151221 The ratio of the gap to the energy spread 
of the ground-state manifold is larger than that of the 



Haldane model. As shown in Fig. 23 the energy gap AE 
remains open and scales to a finite value in the limit of 
N — » oo with N x /N y finite. And the three degenerate 
ground states exhibit spectral flow upon flux insertion. 
The period of 3 fluxes, shown in Fig. [24j indicates the 
system has Hall conductance a xy — 1/3. 

We probe the quasihole excitations by the particle en- 
tanglement spectrum of the ground state. In Fig. 25 
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FIG. 23: Energy gap AE of the Kagome lattice model for 
different system sizes and aspect ratios. In each case, N x = 
3N/N y . 
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FIG. 24: Spectral flow of the low-lying states of the Kagome 
lattice model with N — 10 particles on the N x x N y = 6x5 
lattice upon flux insertion along the x direction. y x counts 
the number of fluxes inserted. The 3-fold ground states flow 
into each other, and do not mix with higher states during flux 
insertion. It takes 3 full fluxes for the 3-fold states to return 
to the original configuration (inset). 



we observe a clear, very large entanglement gap in the 
spectrum, and the counting of the entanglement energy 
levels below the gap again matches in ea ch momentum 
sector the (1, 3)-admissible countin g ^ 15 * 22 !. The width of 
the entanglement gap is A£ = 4.64. This means that 
in the ground state of the Kagome lattice model, the 
inter-particle correlations that obey the (1,3) generalized 
Pauli principle is stronger than any other correlations by 
two orders of magnitude! We thus conclude that the 
ground state is a FCI state with characteristics of the 
FQH Laughlin state. 

Now we briefly address the issue of parameter depen- 
dence. Without loss of generality, we fix t\ = 1 and vary 
Ai in the range (0,V3)- In this region, the single-particle 
spectrum is gapped and the lowest band has unit Chern 
number 32 . A strong correlation between the energy gap 



FIG. 25: Particle entanglement spectrum of the ground state 
of the Kagome lattice model of N = 12 particles on the 
N x x N y = 6x6 lattice, with Nb = 7 particles traced out. The 
number of states below the gray line is 2530 in each momen- 
tum sector, in agreement with the (1, 3)-admissible counting 
rule. The width of the entanglement gap is A£ = 4.64, mean- 
ing that spurious, non-FQH correlations have a probability of 
e -4.64 ^ 0.0097 relative to the universal ones. 



< 




FIG. 26: The correlation between the standard deviation as 
of Berry curvature and the gap AE in the energy spectrum of 
8 particles on the N x x N y = 6x4 lattice. The calculations are 
performed at t\ = 1. The color of each scatter point encodes 
the value of Ai per the color bar to the right. A negative 
correlation between an and AE is visible. 



AE and the Berry curvature fluctuation 2tto~b on Ai is 
clearly visible in Fig. [26] 

As a sanity check, we have also looked at the case where 
one partially fills the second band. We assume the low- 
est band is completely filled and inert, and diagonalize 
the second band directly. In such a situation, there is 
no evidence for Laughlin-like physics at filling 1 + 1/3. 
This is expected since the second band has zero Chern 
number^. 

In Ref. |32l it was shown that by adding properly cho- 
sen NNN hopping terms, the band gap to bandwidth 
ratio can be enhanced by more than an order of mag- 
nitude, and it was argued that this alternative model 
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FIG. 27: Particle entanglement spectrum of the ground state 
of the Kagome lattice model with NNN hoppings of N = 12 
particles on the N x x N y = 6x6 lattice at (\i,\2,t2) = 
(0.28, 0.2, —0.3), with Nb = 7 particles traced out. The num- 
ber of states below the gray line is 2530 in each momentum 
sector, in agreement with the (1,3) counting rule. The width 
of the entanglement gap is A£ = 1.95, much smaller than the 
gap of the model with only NN hoppings shown in Fig. |25| 



could support a more robust Laughlin-like phase when 
interactions are added. Unfortunately, we observe the 
opposite effect: compared with the simple model with 
only NN hoppings, both the energy gap and the entan- 
glement gap are only about half as large for the two cases 
with NNN hoppings studied in Rcf. 321 The reduction 
of A£ from 4.64 to 1.95 means the relative strength of 
the (1,3) Pauli principle exclusion is reduced by a factor 
of ~ 15. The energy and the entanglement spectra for 
the case with the largest band gap to bandwidth ratio 
[at (Ai,A 2) i 2 ) = (0.28,0.2,-0.3)] are shown in Figs. [22 
and [27] We have checked the nearby parameter region 
and no qualitative change is observed. Carefully tuning 
the NNN hoppings could achieve ~ 10% natter Berry 
curvature, but the energy gap and the entanglement gap 
still end up smaller than the simple model with only NN 
hoppings, showing that the flat Berry curvature - FQH 
phase correspondence is to be taken as a general trend, 
rather than a quantitative result. 



B. Half filling 

Among the FQH phases, the Moore-Read (MR) state 
plays a special role. It is one of the best candidates to 
explain the experimentally observed fraction v = 5/2 and 
its excitations obey non-Abelian statistics. The Z n Read- 
Rezayi (RR) states are the generalization of Laughlin 
(n = 1) and MR (n — 2) states. They occur at a filling 
factor v = n/(n + 2). The Z3 RR state is a potential 
candidate to explain the fraction v = 12/5. 

It has been recently shown that a (n + l)-body short- 
range interaction can stabilize the analogue of a Z n RR 
state on the checkerboard lattice^. The (n + l)-body 



interaction used in Ref. [55J mimics the short range inter- 
action for which the Z n RR state is the exact zero-energy 
ground state with highest density. It has also been very 
recently demonstrated^ that the Haldane model of hard- 
core bosons exhibits a MR state at filling v = 1. The 
three-body hard-core repulsion in that case was imple- 
mented by restricting the occupancy of any site to be 
less than three. 

In this section, we focus on the half-filling case and 
thus on the fermionic MR state. Until now, we have 
found no evidence that a two-body NN interaction could 
stabilize a MR-like state in any of the lattice models stud- 
ied. Instead of the two-body NN repulsion, we add to the 
Kagome lattice model a three-body NN repulsion similar 
to the one in Ref. [55J as shown by the gray triangles in 
Fig. [2l] After the gauge transform, the density-density- 
density interaction term reads 



^ E S ™tl +k 3 ^Ub^Lc^sC^b^a 

— k 4 — k 6 

X Vkjkakakiksks, (10) 



N 



where 

V4 lk2 k 3 k 4 k 5 k 6 = 1 + e-^-^i-'Cks-k.)-** (11) 

It is well known that a short-range three-body repul- 
sion stabilizes the Pfaffian MR FQH state^ in a half- 
filled Landau leveP^. We therefore expect an analo- 
gous FCI phase appearing at half filling. The MR FQH 
state has fractional excitations governed by the general- 
ized Pauli principle of having no more than 2 particles 
in 4 consecutive Landau level orbitals and the fermionic 
stat istics of no more than one particle allowed per or- 
bita l 23 * 24 4 This enables a direct extension of the (1,3)- 
admissible FCI counting of having no more than 1 parti- 
cle in 3 consecutive orbitals^to a (2, 4)-admissible count- 
ing^ for the possible MR FCI state. 

We diagonalize the interacting Hamiltonian in the flat- 
tened lowest band at filling 1/2. We show the energy 
spectrum of N = 8, 10, 12 particles on the N x x N y — 
In the three cases, a 6-fold degen- 
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~ x 4 lattice in Fig 
erate ground state Is seen at total momenta {6 x (0,0)}, 
{(0, 0), 2x (0, 1), (0, 2), 2x (0, 3)}, {2x (0, 0), 4x (0, 2)}, re- 
spectively. This agrees exactly with the (2, 4)-admissible 
counting proposed above. As shown in Fig. [29] the en- 
ergy gap AE remains open and scales to a finite value 
in the limit of N — > 00 with N x /N y finite. The six de- 
generate ground states exhibit spectral flow upon flux 
insertion. As shown in Fig. |30[ the period of the spectral 
flow is 2 fluxes, and therefore the ground state has Hall 
conductance a xy = 1/2. 

We probe the quasihole excitations by using the par- 
ticle entanglement spectrum of the ground state. In 
Fig. [3Tj we observe a clear and large gap in the en- 
tanglement spectrum, and the counting of the entangle- 
ment energy levels below the gap again matches in each 
momentum sector the (2, 4)-admissible counting as pre- 
dicted in Ref. [551 The width of the entanglement gap is 
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FIG. 28: Low energy spectrum of the Kagome lattice model 
with three-body interactions of TV = 8 (marked by red 
crosses), TV = 10 (green plus signs), and TV = 12 (blue trian- 
gles) particles on the N x x N y — ~ x 4 lattice, with energies 
shifted by Eo, the lowest energy for each system size. We 
only show the lowest excited level in each momentum sector 
in addition to the 6-fold ground state. 
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FIG. 29: Energy gap AE of the Kagome lattice model with 
three-body interactions for different system sizes and aspect 
ratios. In each case, N x = 2N/N V . 



A£ = 1.66. This indicates that the inter-electron corre- 
lations that obey the (2,4) generalized Pauli principle, 
i.e. the pairing of electrons, is 5 times stronger than any 
other kind of correlations. Therefore, we conclude that 
the ground state of the Kagome lattice model with the 
three-body NN interactions at half filling is indeed a MR 
FQH state. 

Based on the results at filling 1/3 and 1/2, we make 
the conjecture that under desirable conditions, a short- 
range (n + l)-body interaction could stabilize Z„ RR 
parafermion FQH state^ in a Chern band at filling 
n/(n + 2). Further work is needed to confirm the cases 
n > 3. The n — 3 case has recently been reported in 
Ref.CH 

We have also looked for a MR phase in the Haldane and 
two-orbital models using three-body interactions. Unfor- 



FIG. 30: Spectral flow of the low-lying states of the Kagome 
lattice model with three-body interaction of TV = 10 particles 
on the N x x N y = 5x4 lattice upon flux insertion along the 
y direction. 7 a counts the number of fluxes inserted. The 6- 
fold ground states flow into each other, and do not mix with 
higher states during flux insertion. After insertion of every 2 
full fluxes, the 6-fold states return to the original configuration 
(inset). 
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FIG. 31: Particle entanglement spectrum of the ground state 
of the Kagome lattice model with three-body interactions of 
TV = 12 particles on the TV^ x TV y = 6 x 4 lattice, with TV S = 6 
particles traced out. The number of states below the gray line 



is 2910 in momentum sectors with k x — 1,3 and k v 



0,2, 



2912 in sectors (0,0) and (0,2), 2940 in sectors with k x = 5 
and sectors with k x =2,4 and k y — 2, 4, and 2944 in all the 
other sectors, in agreement with the (2, 4)-admissible counting 
rule. The width of the entanglement gap is A£ = 1.66. 



tunately, these two models do not seem to exhibit such 
a phase. This appears to be consistent with the more 
fragile Laughlin-like phase in these models. 



V. RUBY LATTICE MODEL 

The last Chern insulator we have analyzed is on a ruby 
lattice. Hu et alP^ have considered a two-dimensional 
ruby lattice with strong spin-orbit coupling. The sim- 
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FIG. 32: The Ruby lattice model. The six sublattices 1 to 
6 are colored respectively in red, cyan, green, yellow, blue 
and magenta. The lattice translation vectors are bi, b2. The 
complex hopping parameters between NN are t = t r + it« for 
sites having the same parity (black arrows) and ti = t\ r + 
itu for sites having opposite parity (purple arrows), both in 
the direction of the arrows. The real hopping parameter on 
the diagonal of the square (olive dashed lines) is given by t& . 
The density-density repulsion between NN is depicted by gray 
ellipses. 



plified spin-polarized version of this model was shown to 
provide a Chern insulator with an extremely fiat lowest 
band. 



As depicted in Fig. [32j the ruby lattice is spanned by 
the translation vectors bi and b2 and it is made of six 
sublattices, denoted from 1 to 6. After a Fourier trans- 
form and a trivial gauge transform, the single-particle 
Hamiltonian can be cast in Bloch form as 



H = 



EE 

k i,j'=l 



(12) 



Here the lattice momentum k = (k • b 1; k • b 2 ) = (k x , k y ) 
is summed over the first Brillouin zone, and the h(k) 
matrix is given by (the upper triangle can be obtained 
from Hermiticity) 



h(k) = - 





r* 
H 

t 

t 4 (l + e lk * 
t* 





t 

t 4 (l + e- i ^ +fc »)) 
t* 



^4 (^€? kx -\- g*(^ C;c ^ _ ^y)j 







(13) 



We adopt the parameter values suggested in 
the original article^, namely (£t, iir> iii> ^4) — 
(1.2, — 1.2, 2.6, — 1.2)t r . For these parameters, the 
lowest band of the problem is gapped, and has unit 
Chern number^. 

Again, we flatten the Bloch bands using projectors, 
and we add density-density repulsion of unit strength be- 
tween nearest neighbors. There are 12 type of terms as 
depicted by the gray ellipses in Fig. [32] We diagonal- 
ize the interacting Hamiltonian in the flattened lowest 
band at filling 1/3 with up to N = 12 particles, and 
find results quite similar to the previous three models. 
Namely, we find a 3-fold degenerate ground state in the 
momentum sectors predicted by the (1,3) counting prin- 
ciple. The degenerate ground state is separated from 



the excitations by a finite energy gap. Finite-size scal- 
ing indicates that the gap remains open in the thermo- 
dynamic limit. Twisting boundary drives spectral flow 
within the ground-state manifold and the flow has a pe- 
riod of 3 fluxes. Hence, the system has Hall conductance 
&xy = 1/3. We calculate the entanglement spectra for 
various system sizes and find an entanglement gap with 
the (1, 3)-admissible counting in each case. As an exam- 



ple, we show in Fig. 33 the entanglement spectrum of the 
ground state of N = 12 particles. The entanglement gap 
A£ = 5.07 is comparable to, although slightly smaller 
than the entanglement gap A£ = 5.53 of the Kagome 
lattice model with only NN hoppings at the same system 
size. This pronounced sign of exclusion statistics count- 
ing suggests that the system has fractional excitations 
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FIG. 33: Particle entanglement spectrum of the ground 
state of the ruby lattice model of N = 12 particles on the 
N x x N y = 6x6 lattice, with Nb = 8 particles traced 
out. The number of states below each of the three gray lines 
is 741, 1629, 1645 in momentum sectors where both k x and 
k y are even and 728, 1612, 1628 elsewhere. These numbers 
are in perfect accordance with the (1,3)-, (3,7)-, and (3,5)- 
admissible counting, respectively. The width of the (1,3) en- 
tanglement gap is A£ = 5.07. The counting of levels below 
the unmarked gaps do not correspond to any (k, r)-admissible 
rule. 



similar to the Laughlin quasiholes. We do however stress 
that we have not computed the statistics of the quasi- 
holes, as this would require braiding operations likely to 
be plagued by finite-size effects on the lattices we can 
reach by computers. We conclude that the ground state 
of the ruby lattice model with two-body NN repulsions 
is a FQH Laughlin state. 

We have also checked the effect of three-body NN re- 
pulsions at half filling. There are a total of 14 terms, 
corresponding to 2 equilateral triangles and 12 right tri- 
angles embedded in rectangles. Similar to the case of the 



Kagome lattice model detailed in Section [IV B[ we find a 
robust MR state here. The state has the hallmark 6-fold 
degenerate gapped ground state at the correct momenta. 
We also find a large entanglement gap corresponding to 
fractional excitations governed by the (2,4) generalized 
Pauli principle. 



VI. STRUCTURES IN ENTANGLEMENT 
SPECTRUM 

As discussed in previous sections, the entanglement 
spectrum of the ground state of a short-range (m + 1)- 
body interaction has a gap corresponding to the (m, m + 
2)-admissiblc counting, for m = 1,2. This gap measures 
the prominence of m-particle clustering in the ground 
state. The (m + l)-body interactions could possibly gen- 
erate clusters of other sizes as well. For example, in the 
FQH effect in the continuum, MR states can be obtained 
with just two-body potentials. In that case, there could 



be additional gaps in the entanglement spectrum, and 
the counting of levels below these gaps could conform to 
other generalized Pauli principles. 

The appearance of a (n, n + r)-admissible counting re- 
veals the existence of a specific clustering pattern in the 
ground state. This pattern is also present in the model 
fermionic FQH state at filling v = n/(n + r) described 
by the (n, r) Jack p olyno mials multiplied by a Van- 
dermonde determinant 23 * 24 !. For example, a (n,n + 2)- 
admissible counting corresponds to the Z„ Read-Rezayi 
FQH state^, and a (2, 5)-admissible counting corre- 
sponds to the 'Gaffnian' FQH wave function^. The 
model wave functions of such FQH states have character- 
istic zeros when a cluster of n + 1 particles forms even af- 
ter removing the "trivial" zeros provided by the fermionic 
statistics. This reflects specific (n + l)-body correlations 
in the system. Therefore, the presence of a gap with a 
(n, n + r) admissible counting in the FCI entanglement 
spectrum signals the presence of clustering correlations 
similar to a specific model FQH state, and implies the 
presence of stable (n + l)-body correlations in the sys- 
tem. 

We check comprehensively all the four models stud- 
ied i n this paper as well as the checkerboard lattice 
modefMEII, a t various system sizes up to N — 12 par- 
ticles. In most cases, we find extra entanglement gaps 
other than the one above the states of (m, m + 2)- 
admissible counting in the ground state of an (m + 1)- 
body interaction, and the counting of entanglement en- 
ergy levels below some of these gaps matches with an 
(n, n + r)-admissible counting in each momentum sector. 
For all the systems we look at, we usually find at least one 
such extra gap that matches perfectly with a particular 
(n, n + r)-admissible rule. For example, in the Kagome 
lattice model with two-body interaction, we find two ex- 
tra entanglement gaps in the ground state (Fig. 34), with 
the counting of levels below them given exactly (in each 
momentum sector) by the (2,4)- and (2, 5)-admissible 
counting and the folding based on the FQH-FCI map- 
ping^!. The former case corresponds to the MR state 
while the latter corresponds to the Gaffnian wave func- 
tion. We find the appearance of the counting of non- 
Abelian statistics in the Laughlin-like ground state of 
two-body interactions quite intriguing. 

We do not have a complete understanding of the full 
gap structure at the moment. In the earlier example of 
the ruby lattice model shown in Fig. |33[ we find a few 
extra entanglement gaps. Two of them have counting 
below the gap given exactly (in each momentum sector) 
by (3, 5)- and (3, 7)-admissible rules^l respectively, the 
former of which corresponds to a Z3 Read-Rezayi state. 
There are other entanglement gaps in the spectrum that 
cannot be explained by any (n, n + r)-admissible rule. 

All of such extra (n, n + r) entanglement gaps that 
we observe have n + 1 = Na, where Na is the number of 
particles left in the system when evaluating the entangle- 
ment spectrum. Tracing out part of the system enhances 
the correlation between the remaining Na particles. One 
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FIG. 34: Particle entanglement spectrum of the ground state 
of the Kagome lattice model of N = 12 particles on the 
N x x N y = 6x6 lattice, with Nb = 9 particles traced out. 
The number of states below each of the three gray lines is 
138, 195, 198 in momentum sectors where k x ,k y £ {0, 3} and 
135, 192, 195 elsewhere. These numbers are in perfect accor- 
dance with the (1,3)-, (2,5)-, and (2, 4)-admissible counting, 
respectively. 



can imagine that such enhancement is more significant 
for the correlation between the Na particles as a whole 
than the correlation within a fraction of the Na particles. 
In most cases, we observe such (Na — 1, Na — I+t) entan- 
glement gaps for almost every value of Na.- Our examples 
show that a ground-state wave function at v — 1/3 can 
give us information about possible non-Abelian correla- 
tions. 

This suggests a feature of the FCI: the (n + l)-body 
interacting Hamiltonian could not only generate strong 
(n + l)-body correlations, but also higher-body corre- 
lations. This gives hope for the possible realization of 
other members of the Read-Rezayi series with two-body 
potentials. It is very intriguing that entanglement in the 
ground state of an (n + l)-body interacting Hamiltonian 
at filling v = n/(n + 2) contains clear clues that some 
other FQH states could possibly be stabilized by some 
other interaction at some other filling. 



VII. DISCUSSION AND CONCLUSION 

In this paper, we have studied four examples of flat- 
band Chern insulators. If we include the previous arti- 
cles^HmH] on the checkerboard lattice, we now have five 
cases which can be used to compare the universal physics. 

First, we have showed that the existence of FCI does 
not require a special lattice. Under proper conditions, 
FCI phases could (at least) emerge from lattice models 
with 2, 3, or 6 sites per unit cell. There is no restriction 
on the type of Bravais lattice either. 

Second, we find concrete evidence that the stability of 
the FCI phase is not guaranteed by flat band dispersion. 
A large enough gap to bandwidth ratio is necessary for 



the existence of an insulating state at a fractional fill- 
ing. Our calculations performed in the flat-band limit do 
not capture this effect. However, our results reveal that 
the single-particle eigenstates play a critical role in the 
formation of a FCI state as well; a badly chosen param- 
eter set could destroy the topological phase even when 
the band dispersion is perfectly flat. The example of the 
Kagome model is especially instructive. While the addi- 
tion of NNN hopping terms could make the band flatter, 
it actually deteriorates the FCI phase. Thus the search 
of a realistic FCI that would include the effect of the 
band dispersion relation, should not focus only on the 
flat-band criterion. 

Third, the development of a FCI state is often depen- 
dent on suppressed fluctuations in the Berry curvature. 
This agrees qualitatively with the picture derived from 
the algebra of projected density operators. It should be 
noted, however, that even the suppressed curvature fluc- 
tuations could actually be considerably strong. 

Fourth, out of the five models that we have checked, 
the Kagome lattice model with only NN hoppings has the 
most robust FCI phase. The ruby lattice model comes 
close, but has three more tunable parameters. 

In conclusion, we have established the existence of a 
FCI phase in four distinct flat-band lattice models at fill- 
ing v = 1/3 with short-range repulsive density-density 
interactions. The FCI phase is identified by an incom- 
pressible ground state with Hall conductance a xy = 1/3 
and excitations obeying the same counting as Laughlin 
quasiholes. We have observed that in some cases, band 
structures that favor the emergence of the FCI phase 
have less fluctuating Berry curvature. In the presence 
of short-range three-body repulsive interactions, we have 
also found another FCI phase reminiscent of the Pfaffian 
Moore-Read FQH state at half filling of the Kagome lat- 
tice model and the ruby lattice model, in addition to the 
previously known example of the Moore-Read phase in 
the checkerboard lattice model. This more exotic phase 
has an incompressible ground state with Hall conduc- 
tance o~ xy = 1/2 and excitations exhibiting the count- 
ing of Moore-Read quasiholes. On the technical side, we 
have demonstrated from various angles the power of the 
entanglement spectrum as a sensitive and reliable probe 
of topological order, when no model wave functions are 
available. We have also discussed interesting structures 
in the entanglement spectrum and their implication for 
the possible stable existence of other FQH states at zero 
magnetic field. Future directions for this line of research 
include strongly interacting Chern insulators with Chern 
number higher than one, such as the dice modeP^. 
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